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It is shown analytically that an external tidal gravitational field increases the secular stability of 
a fully general relativistic, rigidly rotating neutron star that is near marginal stability, protecting 
it against gravitational collapse. This stabilization is shown to result from the simple fact that 
the energy 5M(Q,R) required to raise a tide on such a star, divided by the square of the tide’s 
quadrupole moment Q , is a decreasing function of the star’s radius R, (d/dR)[5M(Q, R)/Q 2 ] < 0 
(where, as R changes, the star’s structure is changed in accord with the star’s fundamental mode of 
radial oscillation). If (d/dR)[6M(Q, R)/Q 2 ] were positive, the tidal coupling would destabilize the 
star. As an application, a rigidly rotating, marginally secularly stable neutron star in an inspiraling 
binary system will be protected against secular collapse, and against dynamical collapse, by tidal 
interaction with its companion. The “local-asymptotic-rest-frame” tools used in the analysis are 
somewhat unusual and may be powerful in other studies of neutron stars and black holes interacting 
with an external environment. As a byproduct of the analysis, in an appendix the influence of tidal 
interactions on mass-energy conservation is elucidated. 

04.40.Dg, 04.30.Db, 97.60.Jd, 97.80.-d 


1 


I. INTRODUCTION AND SUMMARY 

Wilson, Mathews, and Maronetti Q have carried out 
fully relativistic numerical simulations of the radiation- 
reaction-induced inspiral of a binary neutron star system. 
To make their computations tractable, they employed 
several approximations of ill-understood accuracy. The 
stars in their simulations were identical and were near 
the maximum allowed mass for an isolated neutron star. 
Correspondingly, when the stars were far apart in their 
orbit, each was stable against gravitational collapse. Sur¬ 
prisingly, as the stars spiraled inward, the simulations in¬ 
dicated that their gravitational interaction destabilized 
them, triggering them to collapse before their inspiral 
ended. The magnitude of the destabilization and mathe¬ 
matical arguments to explain it |jj suggest that it should 
show up in the first post-Newtonian approximation to 
general relativity. 

Several researchers have argued that this surprising 
destabilization is wrong: Lai Q has shown that tidal 
interactions between two nearly Newtonian stars will 
tend to stabilize them against gravitational collapse, not 
destabilize them. Lai’s stabilization effect is formally 
of Newtonian origin, but because of the compactness 
of neutron stars, its magnitude is of much higher post- 
Newtonian order. Wiseman || has elucidated Lai’s con¬ 
clusion by showing that at first post-Newtonian magni¬ 
tude, the stars’ gravitational interactions do not alter 
their individual central densities; and Brady and Hughes 
IQ have shown the same at first order in the mass ratio 
M2/Mi when the two stars are fully relativistic and one is 
much less massive than the other, M2 <C M\. Baumgarte 
et. al. P| have carried out numerical simulations of the 
fully relativistic equilibrium states of a binary neutron 
star system in synchronous, circular orbit — simulations 
analogous to those of Wilson, Mathews, and Maronetti. 
Not only do these simulations show no sign of interaction- 
induced collapse; when combined with “turning-point” 
criteria for secular stability, they actually reveal a stabi¬ 
lization of the stars. 

On the basis of these analyses, it seems likely that the 
destabilization seen by Wilson, Mathews and Maronetti 
does not occur in reality, but instead is an artifact of 
poor-accuracy approximations or is due to some error in 
their computations. 

In this paper we present another analysis that reveals 
stabilization, not destabilization. Our justification for 
yet another paper on this subject is two-fold: First, our 
analysis has broader validity than previous ones—it is 
fully relativistic, not post-Newtonian; and unlike the two 
previous fully relativistic analyses of stability it is 
fully analytic and not based on numerical simulations, it 
permits the stars to rotate with arbitrary angular veloc¬ 
ity (though with spins aligned with the orbital angular 
momentum), and it allows an arbitrary mass ratio. Sec¬ 
ond, our analysis employs an unusual approach, which 
may be useful for other problems in fully relativistic bi¬ 


nary evolution: it is formulated in the local asymptotic 
rest frame of one of the two stars and employs energy 
and angular momentum arguments that relate to that 
star alone and not to the binary system as a whole. Al¬ 
though this approach is unusual within general relativity, 
it is well known in Newtonian and post-Newtonian the¬ 
ory. It, in fact, is a relativistic generalization of Lai’s j^| 
post-Newtonian proof of stabilization. 

This paper is organized as follows: In Sec. II we treat 
an idealized problem that illustrates our method: the 
stabilization of a non-spinning neutron star placed in a 
non-rotating external tidal gravitational field. Then in 
Sec. Ill we generalize to a spinning star and rotating 
tidal field, and as an application we deduce the secular 
stabilization of a spinning neutron star in an inspiraling 
binary. In Sec. IV we argue from this secular stabiliza¬ 
tion result that, if an inspiraling binary’s neutron stars 
are secularly stable at large orbital radii, then they can¬ 
not be dynamically destabilized during the inspiral; and 
we make some concluding remarks. In an appendix we 
elucidate the influence of tidal interactions on energy con¬ 
servation, in the Newtonian limit. The relativistic version 
of this issue is a central aspect of the proof of stability 
given in the body of the paper. 

Throughout we use geometrized units in which the 
speed of light and Newton’s gravitation constant are 
unity. 

II. STATIC STAR AND STATIC TIDAL FIELD 
A. Momentarily Static, Spherical Star 

Consider a family of nonrotating, spherical neutron 
stars characterized by a one-parameter equation of state 
P = P(p), where P is pressure and p is density of to¬ 
tal mass-energy. We shall refer to these stars as equi¬ 
librium configurations and shall denote their masses (as 
measured by the Keplerian orbits of distant planets) by 
M e , their total number of baryons by N e , and their radii 
(defined as circumference/ 27 r) by R e . The mass-radius 
relation M e (R e ) of these equilibrium configurations has 
the qualitative form shown in Fig. |l], and N e (R e ) has a 
similar shape; cf. Refs. |j|| 7 ]. 

As is well known, the equilibrium configuration of max¬ 
imum mass (solid dot in Fig. |l| “ critical configuration ”; 
mass M 0 and radius R 0 ) is secularly marginally sta¬ 
ble: It possesses a zero-frequency mode of expansion (or 
contraction) that takes it horizontally in the figure to 
another equilibrium configuration with the same mass. 
Equilibria on the larger-radius side of the critical con¬ 
figuration are secularly stable; those on the lower-radius 
side are secularly unstable [^] . 
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FIG. 1. The Mass-Radius curve for static, spherical equi¬ 
librium configurations (neutron stars) with some equation of 
state -P(p). 

By secular is meant a mode of stellar deforma¬ 
tion which is slow enough for pycnonuclear reactions 
(pressure-induced nuclear reactions) to keep its matter 
always at the endpoint of nuclear evolution, so the pres¬ 
sure and density changes experienced by the stellar mat¬ 
ter follow the same equation of state P(p) as charac¬ 
terizes the equilibrium configuration itself. The slowest 
of the pycnonuclear reactions are “modified URCA reac¬ 
tions” (essentially (3 and inverse-/? decays), driven as the 
star deforms by the rising or falling Fermi energies of the 
star’s electrons, protons, and neutrons; they can require 
timescales of minutes or longer to equilibrate near and 
below nuclear densities, so in principle the stellar defor¬ 
mations can be secular only on timescales longer than 
this. (In practice, these slow reactions have only a weak 
effect on the equation of state near and above nuclear 
densities, so their slowness is often ignored for near crit¬ 
ical neutron stars.) 

Faster ( dynamical) motions, in which some of the py¬ 
cnonuclear reactions do not go to completion, will be 
characterized by a stiffer equation of state (higher adia¬ 
batic index) than secular motions and thus will be more 
stable. Correspondingly, all equilibria on the large-i? e 
branch, being secularly stable, must also be dynamically 
stable, as must be the critical configuration itself. For an 
ancient, further discussion of these issues but in different 
language, see Ref. || 

In the next subsection we shall study the influence of 
an external tidal field on the secular stability of configu¬ 
rations that are nearly critical (nearly at the maximum 
of the mass-radius curve). As an aid in that study, it will 
be useful to consider momentarily static, spherical stars 
that are deformed slightly away from equilibrium. For 
a star containing N baryons (with N very nearly equal 
to N 0 ), we obtain such a non-equilibrium configuration 
as follows: Begin with the equilibrium configuration that 
has N e = N and radius R e (N e ). Expand it or contract it 
slightly to the desired new radius R. In this deformation, 
displace each fluid element by an amount proportional to 
the equilibrium configuration’s fundamental-mode eigen¬ 
function £(x) of secular vibration, and then remove all ki¬ 


netic energy of deformation. The resulting deformed con¬ 
figuration will have a mass M (N, R) that differs slightly 
from the equilibrium mass M e (N e = N). 

In practice, the mass of the deformed star can be com¬ 
puted using not the star’s true eigenfunction, but rather 
the eigenfunction £„(:?) of the zero-frequency mode of the 
critical configuration. This is because the mass is an ex¬ 
tremum with respect to deformations of the star j(j] ; the 
two eigenfunctions £ and £ 0 differ by fractional amounts 
of order n = (N — N 0 )/N 0 , so the masses of the config¬ 
urations obtained by deformations to radius R via the 
£ motion and the £ 0 motion will differ by a fractional 
amount of order n 2 - which is never of interest in this 
paper. 

The configurations obtained by the above construction 
are characterized by two parameters, (N, R ); and their 
masses are functions of these parameters, M(N,R). 

For ease of analysis, we shall now convert to dimen¬ 
sionless variables: 


to = 


M - M 0 
Mo 


n 


V Y 
N 0 


r = 


R — R 0 


R 


O 


(1) 


These variables characterize a configuration’s fractional 
deviations from the critical configuration. 



FIG. 2. Solid curve: the dimensionless mass-radius relation 
for the equilibrium configurations of Fig. |l[ Dashed curves: 
the dimensionless mass-radius relations m(n, r) for configu¬ 
rations of fixed baryon number n that are obtained from an 
equilibrium configuration via deformation along the funda¬ 
mental eigenfunction of radial secular motion. 


Figure^ shows the dimensionless mass-radius relation 
m(n, r) for configurations with fixed baryon number n 
(dashed curves), along with the equilibrium configura¬ 
tions (solid curve). Because the equilibria on the positive- 
r branch are stable against secular deformations (with 
eigenfunction £), they lie at minima of the dashed curves; 
because those on the negative-r branch are unstable, they 
lie at maxima of the dashed curves. 

These dashed curves, when expressed as a power series, 
have the following form: 

to- = ar 3 + (bo + b\r)n , (2) 

where higher-order terms are of no importance in this 
paper, and where a, bo , and b± are all positive. This 
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form is dictated by the demand that for n = 0 the dashed 
curve be flat at r = 0, and that for n < 0 the positive- 
r branch have a minimum (stable equilibrium) and the 
negative-r branch have a maximum (unstable), and that 
for n > 0 there be no equilibria at all (no extrema of the 
dashed curves). 

The equilibrium configurations are located at the ex¬ 
trema of these mass-radius curves, i.e. at locations where 
( dm/dr) n = 0, which yields the following equilibrium 
relations: 


—3a 


n e = 


m _ — 3abo o _ 2 ar 3 


71l e — boTl e T 2 


bi e ’ h 

(foi/ 3) 3 / 2 


I 1 / 2 


(-».) 3/2 


( 3 ) 


where the upper sign is for the right branch (stable stars; 
lower mass at fixed n e ) and the lower sign, for the left 
branch (unstable stars; higher mass at fixed n e ). 

The coefficients a, bo, and b i are determined as fol¬ 
lows in terms of the equilibrium configurations: As is 
well known—cf. Eq. (28) of Ref. || —the mass-energy 
required to create one baryon and inject it into an ar¬ 
bitrary location in an equilibrium configuration, in lo¬ 
cal thermodynamic equilibrium with the matter there, is 
dM e /dN e = fj.By/1 — 2 M e /R e , where /is is the rest mass 
of one baryon at the star’s surface (1/56 the mass of an 
56 Fe nucleus if the star’s matter has been “catalyzed to 
the endpoint of thermonuclear evolution” |]]). Evaluat¬ 
ing this “injection energy” for the critical configuration, 
switching to dimensionless variables, and comparing with 
(■ dm e /dn e ) 0 = bo [Eq. (|j)], we see that 

bo = (pbN 0 IM 0 W 1 - 2 M 0 /R 0 ~ 0.8 . (4a) 


Performing the same calculation slightly away from the 
critical configuration, we obtain 


h _ p bNo/Rq 
1 “ y/1 - 2 M 0 /Ro 


(4b) 


The remaining coefficient, a, is determined by the curva¬ 
ture (d 2 M e /dR e 2 ) 0 of the equilibrium mass-radius rela¬ 
tion at its critical point; cf. Eq. (j|[): 

- R ° f -d 2 Me \ 1 , , , 

6(1-2 M 0 /Ro)\ dR 2 e ) 0 ~ ( C) 

In Eqs. 0) and later equations in this paper, the numeri¬ 
cal values have been inferred, with uncertainties typically 
no worse than a factor two, from the equilibrium config¬ 
urations for plausible equations of state {l]. 


B. Static, Tidally Deformed Star 

We now place a near-critical neutron star in a static, 
external tidal gravitational field, which we characterize 


by the space-time-space-time components of its Riemann 
tensor, Rjoko = £jk [)§• This tidal field will deform the 
star, i.e. will gravitationally “polarize” it, giving it a grav¬ 
itational mass quadrupole moment Ijk- The tidal field 
£jk, quadrupole moment Ijk, and total stellar mass M 
(including the deformation energy and the energy of in¬ 
teraction between the deformation and the tidal field) 
all show up as coefficients in a power series expansion of 
the spacetime metric in the star’s local asymptotic rest 
frame @0 For example, in harmonic (deDonder) coor¬ 
dinates that are attached to the star’s center of mass, the 
time-time metric component has the form iffi 


g oo = 


_1 + S £_ 2 £ 
-£jfcrPn fe r 2 + ... , 


M 3 
2 — + 3 


Ijk'nPn k 


+ ... 


( 5 ) 


where the first “...” denotes higher-order terms in 1/r 
and the second “...” denotes higher-order terms in r. 
Here, and only here, r is coordinate radius computed as 
though the spatial coordinates xl were precisely Carte¬ 
sian (elsewhere in this paper r = (R — R 0 )/R 0 is the 
star’s dimensionless radius); and iR = /r is the “unit 
radial vector”. The local asymptotic rest frame, where 
the expansion (|5j) is valid, is the region from the neutron 
star’s surface out to a distance where the external tidal 
field can no longer be regarded as uniform. 

For simplicity, and in accord with the case of a tidal 
field produced by a distant binary companion, we shall 
assume that the tidal field is axisymmetric; and we shall 
choose our (nearly Cartesian) harmonic coordinates so 
its symmetry axis is along the x 3 direction. Then the 
induced quadrupole moment also possesses this symme¬ 
try, and the components of the tidal field and quadrupole 
moment take the form 

£xx — £yy — ~^£zz — £, Ixx — dfyy —T 2Z = Q. 

( 6 ) 


Here the signs are chosen such that £ and Q are both 
positive. 

In the next few paragraphs, culminating with Eq. 
we shall compute the magnitude Q of the quadrupole 
moment that is induced by a given tidal field £. We do 
so by the following thought experiment: Begin with a 
near critical, spherical configuration of mass M, baryon 
number N, and radius R , and turn on the tidal field £ 


1 Reference j^] is the principal conceptual and mathemati¬ 
cal foundation for this paper’s analysis. The physical concepts 
that underlie Ref. ]TcJ and this paper, including the validity of 
the equivalence principle for “extended” self-gravitating bod¬ 
ies such as neutron stars, date back to John Wheeler’s discus¬ 
sion of equations of motion in Sec. 20.6 of MTW j^] and to 
references therein. 
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without letting the star deform. Then allow the star to 
deform of its own accord, in response to the fixed tidal 
field £. The deformation will develop on the timescale of 
the star’s /-mode quadrupolar oscillations |] ~ 0.3msec, 
which is far faster than the star’s radial motions; thus, 
its angle-averaged radius R will remain essentially un¬ 
changed during the deformation. As the deformation pro¬ 
ceeds, i.e. as Q grows, the tidal field does work W(Q,£) 
on the star, i.e. it increases the star’s total mass-energy 
(excluding quadrupole-tidal interaction energy) by that 
amount. Some portion 6M(Q,R ) of the work W goes 
into producing the stellar deformation (pushing mass up 
in the polar regions and down in the equatorial regions); 
the rest goes into kinetic energy of quadrupolar vibra¬ 
tions. Suppose that we extract the kinetic energy; then 
the star will settle down into the deformed configuration 
that has released the most kinetic energy, i.e. the config¬ 
uration that minimizes the potential energy function 

V(N, R, Q, £) = M(N, R) + SM{Q , R) - W(Q, £) . (7) 

The quadrupolar deformation energy 5M(Q, R) can be 
deduced in cgs units by dimensional arguments. It ob¬ 
viously must be quadratic in Q, proportional to New¬ 
ton’s gravitation constant G and independent of the 
speed of light c, which means it must have the form 
SM = f3GQ 2 /R 5 where R is the configuration’s radius 
(the only lengthscale available other than those, of order 
R , that characterize the star’s internal structure). Here (3 
is a dimensionless coefficient of order unity that depends 
on dimensionless aspects of the star’s internal structure. 
In geometrized units (G = c = 1), this mass-energy of 
deformation is 

SM(Q,R) = -^Q 2 . (8) 

This is the same expression as one obtains in Newtonian 
theory. It can be understood, in Newtonian language, as 
the gravitational energy ~ (A M) 2 /R of an excess mass 
AM ~ Q/R 2 moved into the star’s polar region a little 
higher than radius R, and a corresponding mass deficit 
—AM in the star’s equatorial region a little lower than 
radius R. As an example, for a uniform-density Newto¬ 
nian star ji(|, p = 3. 

The work W(Q, £), which the fixed tidal field £ does 
on the star as its quadrupole moment grows from zero 
to the value Q, can be deduced by examining the flow 
of energy through the star’s local asymptotic rest frame. 
That energy flow can be described by the Landau-Lifshitz 
pseudotensor 0 for the metric (||) with £ fixed and Q 
time varying. By integrating the pseudotensor over a 


sphere S in the star’s local asymptotic rest frame, we 
obtain the rate of change of the star’s total mass-energy 
[the quantity M appearing in the metric (!)]: 

= - J s (-9)t J LL n j d2A ; ( 9 ) 

cf. Secs. 20.2 and 20.3 of Ref. Jl8[ . This dM/dt consists 
of two parts: the rate ( d/dt)WXQ ,£ ) that work is done 
by the external field on the star, and the rate of change 
dEi n t/dt of the interaction energy of the external field 
and the stellar deformation 

dM = dW(Q,£) dBint (1Q) 

dt dt dt 


This split of dM/dt into two parts is elucidated in the 
Appendix and will be analyzed at greater length in a sub¬ 
sequent paper ]l9j ]. The integral <D has been evaluated 
by Zhang |g(J using techniques described in Ref. @ and 
foundations laid in 1211. The result, in a general situation 
where both £ and Q may be changing, isQ 


dM 

dt 


dW(Q,£) dEfr 


dt 


dt 

o ') jj irud ') / 


= 3£ 


2 J dt 10 dt 
dQ 3 d{£Q) 
dt ^ 5 dt 


( 11 ) 


The following argument tells us how much of this mass- 
energy change goes into work and how much into interac¬ 
tion energy: (i) The interaction energy £i n t can depend 
only on the instantaneous stellar deformation and tidal 
field, so dE- lnt /dt must always be a perfect differential. 
By contrast, the rate dW/dt that work is done need not 
be a perfect differential, (ii) In a physical situation (not 
ours) where the tidal field is changing while the stellar 
configuration is constant, to first-order in the tidal per¬ 
turbation no work is done on the star, dW/dt = 0. These 
two facts are sufficient to imply that 


dW(Q,£) = _1 = 3£ dQ 

dt 2 u dt dt 


( 12 ) 


dE int 1 d{£ijTi) = 3 d(£Q) 

dt 10 dt 5 dt 


(13) 


In our thought experiment, where the tidal field is fixed 
at £ and the quadrupole moment grows from zero to Q, 
the total work done on the star is 


2 The tidal field, having no radial nodes in its tidal force pat¬ 
tern, will primarily excite the node-free /-mode; it will cou¬ 
ple far more weakly to the lower-frequency, longer-timescale 
g-modes 


3 Zhang does not give explicitly the numerical coefficient 
— 1/10, since the perfect time derivative term is not of interest 
for his problem; it can be derived by filling in intermediate 
steps in Zhang’s calculation. 
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W(Q,£) = ~l ij £ ij = 3Q£ . (14) 

By inserting expressions ([bi[) and (||) into (7) and min¬ 
imizing the resulting potential energy V ( N , R, Q, £) with 
respect to Q at fixed N, R , £, we deduce the equilibrium 
value of the quadrupole moment: 



(15) 


Correspondingly, the potential energy of the quadrupolar 
equilibrated configuration is 

Q /? 5 

V(N, R , £) = M (TV, R) - —£ 2 . (16) 


Because the quadrupolar contribution — (9R 5 /4f3)£ 2 to 
this potential energy depends on the star’s radius, the 
tidal coupling influences the star’s radial, secular mo¬ 
tions. To deduce that influence, we convert the potential 
energy to dimensionless units, 


)(n,r,e) = 


V(N, R, £) — M a 


£ 


£ = 


M 0 ’ " M 0 /R 0 3 

and combine with Eqs. © and !) to obtain 
v(n, r, e) = ar 3 + (b 0 + b\r)n — (c 0 + cir)e 2 . 


(17) 


(18) 


Here e is the external tidal field measured in units of the 
tidal field produced by the critical configuration near its 
own surface, and 


co = 


ci = 


1 

Mo 

Rp 



9R a 5 9 M 0 

- 4 Wo 

' d_ /9 R 5 Y 
dR 


0.3 , 


5co 


1.5, 


(19) 


where the factor 5 comes from differentiating the R 5 . The 
dependence of /3 on the star’s internal structure may give 
rise to a slight variation of (3 with R, which may slightly 
change the relationship ci = 5co; hence the approximate 
equality in Eq. ©• 

As the near critical star changes its radius R slowly 
and secularly at fixed tidal field £, its quadrupole mo¬ 
ment continually equilibrates, with an accompanying 
tidal feeding of energy W into and out of the star. 
These radial motions conserve the star’s relativistic ki¬ 
netic energy plus its potential energy V(N,R,£) [which 
includes a correction for the flow of W ; cf. Eq. ©]■ 
Correspondingly, the star’s equilibria are the extrema of 
V(N, R, £)—or equivalently of its dimensionless version 
v(n, r, e) —and these equilibria are stable if they minimize 
v and unstable if they maximize it. 

By differentiating Eq. (|l8|) with respect to r at fixed n 
and e, we obtain for the equilibrium configurations 


The configurations on the + branch ( R > R a ) are stable, 
and those on the — branch ( R < R 0 ) are unstable. 

Notice that, in response to the tidal field, each stable 
equilibrium configuration increases its radius (r e > 0); 
and correspondingly, since the fundamental radial-mode 
eigenfunction £(x) that describes the radial shape of its 
spherical deformation has no nodes, its central density 
goes down. Not surprisingly, this leads to a secular sta¬ 
bilization of the star: whereas before the tidal field was 
turned on, the maximum number of baryons that the star 
could support without collapsing was N 0l afterward the 
maximum baryon number has increased by a fractional 
amount 

n e max = ^£ 2 ~ 2.5 . (21) 

The tidally induced increase of equilibrium radius ( |2(i| ) 
and increase of maximum baryon number ( pl| ) are the 
same in sign and order of magnitude as Lai has deduced 
previously using post-Newtonian arguments Q. 

The inability of fully relativistic gravity to produce 
a secular instability can be traced to the robustly pos¬ 
itive sign of the coefficient c±, which in turn follows 
from the fact that the star’s deformation energy at fixed 
quadrupole moment, SM(Q, R)/Q 2 = (3/R 5 is a decreas¬ 
ing function of radius. It is very hard to imagine any 
neutron star for which this would fail to be true. 


III. ROTATING STAR AND ROTATING TIDAL 
FIELD 


We now turn attention to a rigidly rotating neutron 
star interacting with a rotating external tidal field, such 
as those which occur in binary neutron star systems. In 
our analysis, the rotation is with respect to the star’s 
local asymptotic rest frame, which itself will generally 
rotate with respect to “infinity” due to dragging of iner¬ 
tial frames by the binary’s total angular momentum. 

For simplicity, we shall require that the star’s spin axis 
and the tidal field’s rotation axis coincide, as will be the 
case in a binary if the two stars’ spin axes are perpendic¬ 
ular to the orbital plane. This requirement protects the 
tidal field and the spin axis from processing. 

The tidal stabilization analysis of Sec. || can be gen¬ 
eralized to such a rotating system with little change. 
We shall sketch that generalization, beginning in Sub¬ 
sec. Ill A with the secular stability of the rotating star 
in the absence of the tidal field, followed in Subsec. Ill B 
with the star-tide interaction and tidal stabilization. 


A. Rotating Star Without Tidal Field 


r e = ±\J{l/3a)(-bin e + ci£ 2 ) 

. 0 ^ near criticality, where n e <C —£ 2 . (20) 

MJR 3 


±- 


Consider a family of rigidly rotating neutron stars that 
are characterized by the equation of state P{p) and that 
all have the same spin angular momentum J as measured 
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by frame dragging in their local asymptotic rest frames. 
For such stars, because J is fixed once and for all, the 
equilibrium configurations form a one-parameter family 
analogous to that for static stars. We shall character¬ 
ize those equilibria by their masses M e (as measured by 
Keplerian orbits in their local asymptotic rest frames), 
their number of baryons N e , and their equatorial radii 
R e (defined as their equatorial circumferences/27r). 

These equilibria will have a mass-radius relation 
M e (R e ) with the same general shape as that for non¬ 
rotating stars (Fig. [I]), but with a critical mass M a and 
critical radius R 0 that are larger due to centrifugal sta¬ 
bilization and centrifugal forces; cf., e.g., Figs. 1 and 2 of 
Ref. [g. 

As for static stars, so for these spinning stars, the crit¬ 
ical configuration (configuration of maximum mass) is 
marginally stable to a secular expansion or contraction 
that takes it horizontally in Fig. [j] to another equilibrium 
configuration with the same mass M 0 and angular mo¬ 
mentum J; the equilibrium configurations on the larger- 
radius side of critical are secularly stable, while those on 
the smaller-radius side are secularly unstable. 

In these rotating stars, the motions of interest are sec¬ 
ular in two senses: (i) As for static stars, the motions 
must be slow enough for pycnonuclear reactions to go to 
completion so each element of expanding and contracting 
stellar material follows the same equation of state P{p) as 
characterizes the equilibrium stellar structure, (ii) The 
motions must also be slow enough for viscosity to pro¬ 
duce enough coupling between adjacent mass elements 
to keep the star rigidly rotating. This second condition 
is especially severe; so in practice the secular stabiliza¬ 
tion that we shall prove is of interest primaril y be cause 
of its implications for dynamical stability (Sec. 0 

As for static stars, we construct from our one- 
parameter rigidly rotating family of equilibria a two- 
parameter family of configurations that are expanded and 
contracted away from equilibrium. Our construction pro¬ 
cedure is precisely the same as in the static case, with 
the slow, fundamental secular mode of deformation gov¬ 
erning the expansions and contractions. The resulting 
deformed configurations will obviously have mass-radius 
curves, M(N,R) at fixed N (and forever fixed J), with 
the same general shapes as in the static case: the dashed 
curves of Fig. 0. Moreover, when expressed in dimension¬ 
less variables (hi), the m(n,r) relation near the critical 
point for our rotating family will have the same math¬ 
ematical form, m = ar 3 + (bo + b\r)n, as for the static 
case. The coefficients a, bo, and b\ will be affected by the 
rotation and thus will not take on their static forms ([|), 
but in order to produce the equilibria’s known stability 
properties, they will still all be positive. 


B. Tidally Deformed, Rotating Star 


We now place each of our rigidly rotating, near-critical 
stars- all with the same spin angular momentum J—in 
a rotating, external tidal gravitational field. We assume 
that the star and the tidal field rotate about the same 
axis as seen in the star’s local asymptotic rest frame; 
the star’s rotational angular velocity is u > < 27t x 1000/s, 
and if the rotating field is produced by a binary com¬ 
panion, its angular velocity is that of the orbital mo¬ 
tion, n < 2n x 600/s f§. this case, the tidal field 
is exceedingly unlikely to excite the star’s normal modes 
significantly. This is because (i) significant excitation re¬ 
quires resonant coupling; (ii) the circularly polarized /- 
and p-modes have angular velocities of pattern rotation 
cr > 27r x 1500/s (unless the star is rotating close to cen¬ 
trifugal breakup), which is too large to resonate with the 
driving force except under the most extreme of circum¬ 
stances; and (iii) the low-frequency p-modes, which can 
resonate, have only very weak coupling to the tidal field 
p| pf . With this justification, we shall assume that the 
tidal field does not resonantly excite the star’s normal 
modes. 

The rotating tidal field will, however, raise a nonreso¬ 
nant, rotating quadrupolar tide on the neutron star. If 
the star’s material had zero viscosity, the tide would be 
perfectly aligned with the tidal field; i.e., the axes of the 
star’s rotating quadrupole moment Ijk would be identical 
to those of the tidal field £ ? k- This standard Newtonian 
result must be true also in general relativity since the 
relevant physics is nothing but that of simple harmonic 
oscillators: Whenever an undamped oscillator is driven 
by an off-resonance sinusoidal force, the oscillator’s dis¬ 
placement response is precisely in phase with the force. 
As for an oscillator, so for the star, any small viscosity 
will cause a slight phase lag between excitation and re¬ 
sponse: the orientation of the star’s tidal deformation 
and its quadrupole moment will lag that of the tidal field 
by an angle 


(fl — ui) 
(a — w) 2 t* 


( 22 ) 


where t* is the viscous damping time for the star’s 
quadrupolar excitations. For realistic viscosities, this is 
an exceedingly small lag angle HH; so the quadrupole 
moment is very nearly oriented along the tidal field. 

This agreement of orientations prevents the tidal field 
from torquing the star, so the star’s spin angular momen¬ 
tum J is conserved. 

This conservation of J and alignment of tidal field and 
quadrupole moment enable us to carry over the static- 
star analysis of tidal stabilization (Sec. flB) to the ro¬ 
tating case, essentially without change. The tidal field 
and quadrupole moment can be characterized by scalars 


in the manner of Eq. (Oj); raising 
ergy SM(Q,R) given by Eq. (|l 


the tide requires an en- 
; the tide extracts from 


the tidal field an amount of energy W(Q,£) computable 
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by integrating the Landau-Lifshitz pseudotensor over a 
closed 2 -surface in the star’s local asymptotic rest frame 
in the manner of Eq. |). and that integration must give 
expression ( 0 ) aside from a multiplicative factor that is 
close to unity. This multiplicative factor (call it 1 + 7 ), 
which we have not computed, will arise from the syn¬ 
chronous rotation of the tide and star with velocity lot 
at the star’s surface. It cannot differ from unity by more 
than 7 ~ u 2 R 2 < 0.01 (there is no linear term 7 ~ u>R 
because such a term would reverse sign when the star 
reverses direction of rotation.) As in the static case, the 
quadrupole moment adjusts itself so as to minimize the 
potential energy function V(N, R,Q,£) of Eq. (jj]) ; the re¬ 
sulting Q is expression (|lj|) aside from the multiplicative 
factor 1 + 0 (w 2 r 2 ); and the star’s radius changes slowly 
and secularly in a manner governed by the potential en¬ 
ergy function v(n, r, e) = ar 3 + (bo + b\r)n — (co + cir)e 2 
[Eq. (^ 8 [)]. The only difference is in the numerical values 
of the coefficients a, bo, b±, cq, and Ci- 

As in the static case, the key to tidal stabilization is the 
sign of Ci; and as there, so long as the star’s deformation 
energy at fixed quadrupole moment, SM(Q,R)/Q 2 = 
(/ 3/R 5 ) x [14- 0(w 2 !? 2 )] is a decreasing function of radius 
R, the tidal field secularly stabilizes the star, increases its 
radius, and decreases its central density. It is exceedingly 
difficult to imagine a star for which (when w < 2-7T x 600/s 
as it must be in an inspiraling binary |23| ) this would not 
be the case. 


IV. CONCLUSIONS 


ativistic phase of binary inspiral. Such a star must have 
lived for millions of years to reach this inspiral phase, so 
it must be secularly stable as well as dynamically stable 
before the tidal field of its companion begins to affect it 
significantly. The growing tidal field, as we have seen, 
must increase the star’s secular stability. Since secular 
stability implies dynamical stability, the tidally deformed 
star must also remain dynamically stable. 

This conclusion is not at all obvious from the equations 
of the post-Newtonian approximation, as they are usu¬ 
ally written: in a reference frame that is asymptotically 
inertial at infinity. One can identify in those equations 
terms that appear able to produce tidal instability at 
first post-Newtonian order [0. This fact illustrates the 
superficially misleading character of the post-Newtonian 
equations: Apparently magical cancellations || enforce 
the strong equivalence principle, which is fundamentally 
at the heart of tidal stabilization. 

By contrast, the mathematical tools used in this paper, 
being based on the local asymptotic rest frame of the 
star whose stability interests us, are closely linked to the 
strong equivalence principle and lead rather directly to 
the tidal stabilization. These tools are not widely used in 
relativistic astrophysics. They are worth trying whenever 
one is interested in the behavior of a semi-isolated portion 
of a larger relativistic system—e.g., a neutron star or 
black hole interacting with other bodies [^ 0 |. 

As this paper was being completed, I learned of a 
similar, local-asymptotic-rest-frame analysis by Flana¬ 
gan | 2 Jj], which, however, focuses on equations of motion 
rather than energy considerations. 


In this paper we have proved that tidal fields tend 
to stabilize a star against secular gravitational collapse. 
However, for spinning stars in binaries, secular stability 
is irrelevant because radiation reaction drives the binary 
through the regime of relativistic gravity far too quickly 
for viscosity to keep the rotation rigid in pulsational mo¬ 
tions. The relevant issue in this case is dynamical stabil¬ 
ity- 

As is well known (e.g., Ref. (g4|), secular stability im¬ 
plies dynamical stability. This is so for two reasons. 
First, pycnonuclear reactions go to completion in secular 
motions but not in dynamical motions, and as a result the 
secular motions are characterized by a softer equation of 
state than the dynamical motions and thus are less sta¬ 
ble. Second, dynamical motions of an initially rigidly ro¬ 
tating star produce differential rotation; and the viscous 
coupling that converts that differential rotation into the 
rigid rotation of a secular motion will necessarily extract 
energy from the rotation. This means that, beginning 
with the same equilibrium configuration, the dynamical 
motion must lead to a configuration of greater potential 
energy V than the secular motion, which in turn means 
that the equilibrium configuration is more stable against 
the dynamical motion than against the secular motion. 

Consider a rigidly rotating neutron star in the final, rel- 
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APPENDIX 


In this appendix, we elucidate the energetics of 
quadrupole-tidal coupling, Eqs. (H)-(14), by working out 


its details in the Newtonian limit without rotation. 

We consider a nonroting Newtonian star that is 
spherical, aside from a quadrupolar deformation and 
quadrupole moment Ijk that are aligned with an external 
tidal field £ 3 k = <1 ) e ,jk- Here 4> e is the external Newto¬ 
nian potential, and £jk and Ijk are symmetric and trace 
free. We denote by p, p, v and n the (Newtonian) mass 
density, pressure, velocity and specific internal energy of 
the stellar fluid, and by 


<F = <F 0 + 


(23) 



the Newtonian gravitational field and its split into the 
star’s self field d> 0 and the external field <f> e . 

In Newtonian theory, the total energy density and en¬ 
ergy flux of the stellar fluid plus gravitational field can 
be written as 

e 00 = p(^ 2 + n + $) + ^$ J $ j , (24) 

0 Oj = pv 3 (\v 2 + n + ^ + $) - ■ , (25) 

\ 2 P J 47t 

where we make no distinction between covariant and con- 
travariant spatial indices because our coordinates are as¬ 
sumed Cartesian. It is straightforward to verify that the 
energy conservation law 

e“ f + e 0 jj =o (26) 


is satisfied by virtue of mass conservation + {pv 3 ) j = 
0, the first law of thermodynamics, pdH/dt + pv 3 j = 0, 
the fluid’s Euler equation pdv 3 /dt + p$,j +p,j = 0, and 
Newton’s field equation '1\ J7 = 47rp. Here d/dt = d/dt + 
v 3 d / dx 3 is the fluid’s comoving time derivative. 

[Note: Equations (|24|) and (^5|) entail a specific choice 
of how to localize the system’s gravitational energy. 
Other choices are possible: one can add to © 00 the diver¬ 
gence of r]j = ofixfi y (where a is an arbitrary constant) 
and add to O 0j the time derivative of rp without affect¬ 
ing the law of energy conservation (26) or any of the 
system’s physics. This nonuniqueness of localization of 
gravitational energy also occurs in general relativity; cf. 
Chap. 20 of MTW |Q. We shall discuss its consequences 
below.] 

The external field is purely quadrupolar and source 
free throughout the star and the star’s local asymptotic 
rest frame 


4> e = rV , ® e ,jj = 0, (27) 

and its tidal field £ij evolves with time. The star’s exter¬ 
nal self field is purely monopolar and quadrupolar and 
its source, of course, is the star’s mass distribution 

M 32),-nW . . . 

$ 0 = — -- 3 —outside star, =4np . (28) 

Here r = 5ijX l x 3 is radius and n 3 = x 3 /r is the unit 
radial vector. The star’s mass M is constant in time, but 
its quadrupole moment evolves. 

Now consider the total energy E inside a ball V cen¬ 
tered on the star and larger than the star. This energy 
consists of three parts, the star’s self energy E a , the ex¬ 
ternal field’s energy E e , and the energy E- mt of interac¬ 
tion between the star’s quadrupolar deformation and the 
external tidal field 


E n = 


E r = 


-Bint — 




II -i> 

d 3 x , 


87T 


^0,7^0, 


d 3 a 


pd>e 




r] 3 r — _ f. r.. 

a j, — 


(30) 

(31) 

(32) 


Here the value of E- mt follows from the form and source- 
lessness (27) of the external field, the source equation 
and exterior form of the star’s self field (p8|), and the 
usual volume integral for the quadrupole moment lij = 
f v p(x l x 3 — \r 2 5 l3 )d 3 x. 

The law of local energy conservation (^fi|) guarantees 
that the rate of change of the total energy (E9|) is the 
surface integral of the energy flux over the boundary dV 
of the ball V 


dE 

dt 


0° 3 n 3 r 2 d£l . 


1 dv 


Here dfl is solid an gle. By combining Eqs. 
(p3|) , (p7|), ( [H ) and (^l|), we find that 

dE dE e f 1 ■ 2 

-7r = —rr + / —<$>o,t$o,jn 3 r 2 dQ. 

dt dt J dv 4n 


(33) 

(H), 


+ dt l 10 l3k£jk > 


‘ dt 


(34) 


The first term is the rate of change of the external field 
energy, inside the ball, due to the evolution of the tidal 
field £jk■ The second term is the rate of change of the 
star’s self energy due its own field energy flowing into 
or out of the ball V as its quadrupole moment decreases 
or increases; in the limit as the ball’s radius becomes 
arbitrarily large, this term goes to zero. The third term 
is the rate of change of the interaction energy; cf. Eq. 
([32|). By comparing with Eq. (29), we see that the last 
term must drive changes in the star’s self energy 


dE a 

dt 


1 


/ — 

Idv 47r 


;n 3 r 2 d£l — —Ei-j—Ii- 


dt 


(35) 


In other words, the last term is the rate at which the 
external tidal field does work on the star 


dW(Q,£) 1 d in 

dt 2 y dt 


(36) 


This work rate is in perfect agreement with the fully 
relativistic work rate derived in the text using 

the Landau-Lifshitz pseudotensor to localize the gravi¬ 
tational energy and using deDonder (harmonic) gauge. 

Not so the rate of change of interaction energy. The 
Newtonian analysis gives 


dE int \ = 3 d\£,,£•■') 

dt ) Newton 4 0 dt 


(37) 


E= f Q 00 d 3 x = E o + E e +E int 
Jv 


(29) [Eq. ( jl2| i above]; by contrast, the relativistic analysis 
gives this same expression but with a coefficient —1/10 
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rather than +3/10; see Eq. dD- This disagreement does 
not arise from any fundamental difference between New¬ 
tonian theory and general relativity. Rather, it arises 
because the localization of gravitational energy is non¬ 
unique, both in Newtonian theory and in general rela¬ 
tivity [see the paragraph preceding Eq. @)]; and E- mt is 
that portion of the system’s total gravitational interac¬ 
tion energy which resides in the vicinity of the particular 
neutron star on which we are focusing. For a binary 
system, our two analyses distribute the gravitational en¬ 
ergy differently between the stars’ vicinities and the in¬ 
terstellar region, and they thereby give rise to different 
values of E[ nt . If, in the Newtonian analysis, we were 
to change our localization of the gravitational energy via 
the transformation described in the paragraph preceding 
Eq. (^?j), we would alter Ei nt , while leaving unchanged 
the uniquely defined work W done by the tidal field on 
the star of interest. Similarly, if, in our general relativitis- 
tic analysis, we were to change our energy localization by 
switching from the Landau-Lifshitz pseudotensor to some 
other pseudotensor, or by performing a gauge change on 
the gravitational field, we thereby would alter £) nt but 
leave W unchanged. 

This localization dependence of E- m t and uniqueness of 
W will be analyzed in detail in a subsequent paper []l9|. 
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